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Abstract 

Vibrational spectra of closed shell Lennard-Jones icosahedral and cuboctahedral clusters are 
calculated for shell numbers between 2 and 9. Evolution of the vibrational density of states with 
the cluster shell number is examined and differences between icosahedral and cuboctahedral clusters 
described. This enabled a quantum calculation of quantum ground state energies of the clusters 
in the quasiharmonic approximation and a comparison of the differences between the two types of 
clusters. It is demonstrated that in the quantum treatment, the closed shell icosahedral clusters 
binding energies differ from those of cuboctahedral clusters more than is the case in classical 
treatment. 
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I. INTRODUCTION 



The interest of the condensed matter community in cluster physics has been recently 
revived due to exciting technological possibilities offered by new materials in which clusters 
play the role of basic "building blocks" (see e.g. Ref. jlj]). Theoretical investigations of model 
clusters, i.e. assemblies of particles interacting via bina ry ( mostly Lennard- Jones) potential 



has a long history 
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ll2 , [l^, [l^l . Most studies performed so far 



have concentrated on the search of the most stable configuration of iV particles. Typically, 
total potential energy of a cluster is written as a sum over binary interactions of all pairs of 
particles within the cluster. Then, a search for a set of coordinates (or configuration of the 

n 

cluster) which minimize the potential energy is performed [15j. Depending on the number 
of particles within a cluster, very different cluster shapes can be obtained following this 
procedure M- The approach just sketched can require formidable numerical optimization 
procedures [ll|. However, it is completely classical since it produces an absolute minimum in 
the potential energy, rather than the energy of a ground quantum state of the system which 
must be smaller in absolute value. The difference between the two values is the (quantum) 
zero-point energy of the cluster. 

The clusters investigated in this article are known as closed shell clusters. These 
appear for "magic values" of N, and the specific sequence to be investigated here is 
N = 13,55,147,309,561,923.... The characteristic sequence of "magic" numbers has also 
been observed experimentally (see e.g. Ref. [lij]). Closed shell clusters of interest to the 
present article appear in two configurations, icosahedral and cuboctahedral (see Fig. [TJ). 
These clusters can be thought of as being assembled by adding closed " shells" of atoms to 
a single atom located at the origin, i.e. the center of the cluster. The thus "assembled" 
structure can be characterized with the maximum shell number n (see Fig. The two 
types of clusters have quite different geometries, but their common feature is that they have 
the same number of atoms for the given maximum shell number n |7j. This is a very con- 
venient feature which enables a direct comparison of the various physical properties of such 



clusters 



13j. Cuboctahedral clusters can also be visualized as pieces of a crystal with 



face-centered-cubic (FCC) packing. A FCC crystal can be obtained from cuboctahedral 
clusters in the limit of infinite shell number. The same is not true for icosahedral clusters 
which are therefore "noncrystal" 0. Closed shell icosahedral configuration has been shown 
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to have lower total potential energy for clusters smaller than about 9000 particles Q]. The 
actual number of particles (or the shell number) \^ above which cubqctahedral clusters have 



lower potential energy has been a subject of debate (see e.g. Refs. 

BUI 17(). This is not 

important for the purpose of this article since I am going to consider relatively small clusters 
(up to nine shells or 2869 atoms) for which icosahedral configuration is a classically preferred 
one. A plausible reason for the fact that the closed shell icosahedral (CSIC) configuration 
is classically more stable from the closed shell cuboctahedral (CSFCC) configuration is that 
the arrangement of particles which are on the cluster surface is "tighter" in the CSIC than 
in the CSFCC ordering of the cluster (see Fig. [T]). 

The aim of this article is twofold. First, a detailed microscopic calculation of cluster 
vibrations shall be performed, depending on the closed-shell cluster size. It is of interest 
to see whether the different cluster geometries reflect themselves in the cluster vibrational 
properties. The differences in vibrational properties could be exploited to discriminate 
between the different cluster geometries. Some early studies 

Baa 

have dealt with the 

vibrational frequency spectra of clusters, but these studies were limited to quite small clusters 
and there were no attempts to compare the vibrational spectra pertaining to clusters with 
the different symmetries. The results should also be relevant to the studies of vibrations in 
nanoparticles deposited on substrates |18|. 

The Second aim is to reexamine the cluster stability from the point of quantum mechan- 
ics, i.e. to investigate the differences between the CSIC and CSFCC clusters with respect 
to their quantum ground state energy and temperature dependent vibrational entropy. The 
quantum approach is quite easy in the case in which the cluster dynamics can be adequately 
represented by harmonic vibrations, but it should be kept in mind that such treatment is 
adequate only for sufficiently low temperatures. This is the region of temperatures that is of 
interest to this article. There are several points that come to mind regarding the second aim. 
First, one could assume that the fact that the surfaces of CSFCC clusters are less densely 
packed implies that their characteristic vibrations will be "softer", i.e. of lower frequency. 
Thus, zero-point energy of CSFCC clusters could be expected to be smaller, and the shift of 
the ground state energy from absolute minimum of the potential energy lesser than in CSIC 
clusters. This would suggest that quantum ground state energy may even be larger (in abso- 
lute units) in CSFCC clusters, which would promote them to thermodynamically preferred 
configurations of N atoms at zero temperature, T = K. Second, even if the zero-point 
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energy is not sufficiently different in the two cases, in order to find the thermodynamically 
preferred configuration at finite temperature one should minimize free energy of the system. 
As the cluster free energy depends on the features of the cluster vibrational spectrum and 
temperature, it may happen that at some finite temperature, T c , the free energy of the 
CSFCC configuration becomes smaller than the free energy of CSIC configuration, even if 
its ground state, binding energy was smaller than in the CSIC configuration. Third, as 
vibrational frequencies depend on the mass of particles in the cluster, M, quantum effects 
and thermodynamical considerations must be mass dependent. Classically, the cluster sta- 
bility considerations depend only on the binary interaction potential. Thus, any two clusters 
composed of different isotopes of the same element will always have the same configuration. 
Obviously, this needs not be the case in a quantum treatment and quantum effects are going 
to be larger for clusters composed of lighter particles. This article aims at examining the 
three points mentioned above. The results of the article are not directly applicable to real 
clusters, although the model of the cluster adopted could be used to describe noble gas 
clusters and could serve as a point of departure for setting up more complex interaction 
models There are some reservations, however, since it is known that the form of 

the binary potential employed in calculation can influence the cluster optimal shape and its 
dynamical properties |8l Il3j|. 

The article is organized as follows. In Section m I shall briefly describe the adopted 
theoretical approach (quantum quasiharmonic approximation). Section HTT1 deals with the 
model clusters in which the particles interact via Lennard- Jones binary potentials. Vibra- 
tional spectra of CSIC and CSFCC clusters are calculated for shell numbers between two 
(55 atoms) and nine (2869 atoms). The effects of mass and temperature on the Helmholtz 
free energy of clusters are considered on the example of clusters of Ne, Ar and Xe atoms. 
Section HVl summarizes and concludes the article. 



II. PREREQUISITES FOR THE CALCULATION OF THE VIBRATIONS AND 
VIBRATIONAL FREE ENERGY OF CLUSTERS 

The dynamical behavior of a cluster is described by a set of coordinates {ri,T2, ...,rjy} 
which are treated as time dependent variables. The particles within the cluster are assumed 
to interact via a binary potential, v , which depends only on their relative positions, i.e. 
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v(Ti, Tj) = v(\ri — Tj\). The total potential energy of the cluster, V p is assumed to be given 
by a summation of binary interactions over all the pairs of particles in the cluster, 

V p {{Ti 1 ...,T N }) = y £v(\r i -T j \), (1) 

i>j 

where the dependence of the potential energy on the cluster configuration has been empha- 
sized. The configuration of a cluster which minimizes the potential energy function can be 
denoted by a set of coordinates {if, r®, r° N } which denote mean, static positions of the N 
particles within the cluster. Assuming that = — if; |iij| <C |r° — r°|,Vi, j G {1, ...,N}, 
i.e. that the displacements of particles from their equilibrium positions are small, one can 
expand the total potential energy of the cluster in the Taylor series up to the second order 
as, 

Vi({r 1 ,..,M) = y p ({rS,..,r^}) + iEE^(^^] > ( 2 ) 

where a and /3 denote the Cartesian components (x, y, and z) of the vectors. The first 
derivatives of V p with respect to atom coordinates are assumed to vanish, i.e. the cluster is 
assumed to be in a minimum potential energy configuration. This is the well known harmonic 
approximation and it serves as a starting point for the calculation of the cluster normal modes 
of vibration Q], i.e. a set of linear combinations of {u 1; ...Ujv} variables (or eigenmodes), 
each of which corresponds to a vibration of the system with a single frequency j^. The 
Hamiltonian of the problem when written in terms of normal mode coordinates represents a 
set of independent harmonic oscillators whose both quantum and classical dynamics are well 
known. There are 3iV — 6 such oscillators with characteristic frequencies uj p ,p = 1, 3iV — 
6. Six degrees of freedom that do not represent vibrations are three rotations and three 
translations of the whole system. 

Once a set of eigenfrequencies is calculated, one can proceed to calculate the Helmoltz 
free energy of the cluster, F, which is given by 

F = -k B T\nZ } (3) 

where k B is the Boltzmann constant, T is temperature, and Z is the quantum partition 
function of a system of 37V"— 6 independent oscillators. In terms of the eigenmode frequencies, 
one can write 



(4) 



p=i z p=i 



( hu} p 

exp 1 



knT 



where 

V? = Vp({rS I ...,r§ r }), (5) 

is the minimum of classical potential energy of the cluster (classical ground state energy), and 
H is the reduced Planck constant. The sum of the first two terms in Eq. (j3J) represents the 
quantum ground state energy of the cluster, Eq, calculated in the harmonic approximation. 
At constant temperature, the state which represents thermodynamical equilibrium of the 



system is the one which minimizes the Helmholtz free energy |2l|. Note that even at zero 
temperature, the cluster free energy has a quantum, zero-point energy contribution [second 
term in Eq.(j3J)], in agreement with the ideas put forth in the Introduction. Thus, even at 
zero temperature, the thermodynamical equilibrium state of the cluster need not be the 
same as the state which minimizes the classical potential energy of the cluster This of 
course depends on the nature of particles in the system, their mass in particular as is known 



from the studies of systems of He atoms 



221. 



III. POTENTIAL ENERGY MINIMA, VIBRATIONAL SPECTRA AND ZERO- 
POINT ENERGIES OF LENNARD-JONES CSIC AND CSFCC CLUSTERS 

In order to perform a normal vibrational mode calculation, one first has to find the set of 
coordinates {r®,r9,, ...,r^} which make the functional V p stationary in the 3 N- dimensional 
configurational space. If the normal mode calculation were performed in a nonstationary 
configuration, some of the normal mode frequencies would turn out imaginary. This would 
signify the instability of the cluster structure [lflj ]. 

The atoms were assumed to interact via binary Lennard- Jones 6-12 potentials, 



v(\ri - Tj\) = v(r) = 4e 



a\ 12 (a 



(6) 



and were initially arranged in a configuration which has the symmetry of either cuboc- 
tahedral or icosahedral cluster with the nearest-neighbor atom distances set close to the 
Lennard- Jones range parameter a. The thus obtained initial configuration was then allowed 
to relax to a configuration in which the forces acting on each of the atoms in the cluster 
were smaller than some predefined and arbitrarily small absolute force. To obtain the re- 
laxed configuration, each of the atoms in the cluster was displaced a certain distance in the 
direction of the total force acting on it. This was repeated until the absolute value of the 



force averaged over all cluster atoms dropped below predefined force f c . In each step of 
this iterative procedure, the lengths of vectors the atoms were moved along were reduced 
or enlarged, depending on the vectorial characters of forces in a given and preceding step 
actingon a particular atom. Apparently similar relaxation algorithms were used in Refs. 
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]. The value of f c used in the calculations presented below was 2.0 * 10 -13 e/a, 
except for the n = 2 CSFCC cluster for which f c = 2.0 * 10~ 7 e/a was used. For smaller 
values of f c , the algorithm used relaxes the cluster of n = 2 CSFCC initial symmetry to the 
n = 2 CSIC configuration. Whether this signifies instability of CSFCC n = 2 cluster, as is 

n 

the case for n = 1 CSFCC cluster j2j , or is simply a manifestation of inability of the algo- 
rithm to reach presumably shallow and/or narrow minimum in the configurational space, is 
not clear and is of limited importance to this work - the n = 2 CSFCC cluster was treated 
as stable and the calculation of vibrations has been performed. It suffices to say that the 
potential energy obtained for n = 2 CSFCC cluster in this study agrees perfectly with the 
result of the previous study (see Table HJ). Nevertheless, it is also of interest to note that 
the results of Refs. 0, |^ also indicate the instability of n = 2 CSFCC cluster. 

The classical potential energies of the clusters, V® obtained using the procedure explained 
in the previous paragraph are displayed in Table |l] for closed-shell clusters with maximum 
shell number n. These energies depend only on e, and this was used as a scale for V®. 
The minimum nearest neighbor distance within the cluster, r min = min{|rj — = 
1, N, i 7^ j, is also displayed. The values of V® are equal to those obtained in Ref. p| in 
the first six significant decimal places. This confirms the validity of the relaxation algorithm 
used in this work. The minimal nearest neighbor distance found in clusters is always between 
the central atom and one of the atoms in the first cuboctahedral or icosahedral shell. It is 
of interest to note that the nearest neighbor distance in Lennard- Jones 6-12 FCC crystal is 
1.09017a H- 

When the equilibrium configuration was obtained, the force-constants [second derivatives 
of the binary potential function, see Eq.Q] acting between all pairs of atoms within the 
cluster were calculated, the dynamical vibrational matrix was set up and its diagonalization 
performed in order to obtain eigenmode frequencies and polarization vectors [3]. The 
frequencies obtained in such a way depend on the atom mass, M, and on the Lennard- 
Jones potential parameters, but only through their combination, u = J~ej {Ma 2 ) [2^ . This 
combination of parameters was used as a universal frequency scale. 
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TABLE I: Classical potential energies (V p °) in units of e and minimum nearest neighbor distances 
(r m i n ) in units of a of Lennard- Jones icosahedral and cuboctahedral closed-shell clusters. The 
maximum cluster shell number is denoted by n, and total number of atoms in the cluster by N. 



Icosahedral, CSIC Cuboctahedral, CSFCC 
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N 


t 7-0 r i 


Train \p\ 


t 7-0 r 1 


Train [c] 


9 

Li 


55 


-97Q 948^ 


1 05045 


-968 9765 


1 0Q3Q9 


3 


147 


-876.4612 


1.03548 


-854.3766 


1.09093 


4 


309 


-2007.219 


1.02596 


-1971.561 


1.08998 


5 


561 


-3842.394 


1.01904 


-3792.097 


1.08929 


6 


923 


-6552.723 


1.01361 


-6488.217 


1.08894 


7 


1415 


-10308.89 


1.00914 


-10232.14 


1.08872 


8 


2057 


-15281.55 


1.00535 


-15196.07 


1.08842 


9 


2869 


-21641.35 


1.00205 


-21552.22 


1.08821 



Figures El and El display the vibrational densities of states per atom, p(uS)/N of CSIC 
and CSFCC clusters, respectively, as a function of maximum cluster shell number. The 
vibrational density of states was calculated as 

3N-6 

PH = Yl 8(u-u p ). (7) 
P =i 

For the sake of easier visualization, the ^-functions in the above equation were represented 
by normalized gaussians with a width parameter of 0.02 ujq. The phonon density of states 
per atom of FCC Lennard- Jones 6-12 crystal is displayed in the bottom panels of Figs. 
121 and El for comparison. The phonon density of states for FCC crystal was calculated 
by numerical sampling the of three-dimensional inverse (wave vector) space in 3000000 
randomly distributed points. 

One can immediately note quite different scales on x-axes in Fig. |2]and Fig. El This is due 
to the fact that the vibrational spectra of CSIC clusters contain a high-frequency tail which 
is not present in the CSFCC case. For example, the maximum frequency of n = 3 CSIC 
cluster is 0J max = 37.20u; , while in n = 3 CSFCC cluster it amounts to u) max = 25.96a;o- 
The reason for such a large difference can be found in the data presented in Table [1] From 
the inspection of minimum nearest-neighbor distances, one can conclude that the CSIC 
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clusters are more tightly packed and that some of the neighbors are in a repulsive region of 
their binary interaction potential (the minimum of the Lennard- Jones 6-12 binary potential 



is at r = 1.122(7 



24j). This especially applies to the central cluster atom which is very 



tightly surrounded by the atoms in the first icosahedral shell (see Table U}. The modes that 
dominantly represent relative motion of such atoms have therefore quite high frequencies. To 
substantiate this claim, I have plotted in Fig. H]the eigenvectors (or displacement patterns) 
of some of the characteristic modes pertaining to CSIC clusters. The displacement pattern 
of the highest frequency mode of n = 2 CSIC cluster is represented in panel (b) of Fig. 0] 
It can be seen that in this mode the motion of the central atom dominates the displacement 
pattern. Atoms in the first shell also slightly move, so that the total linear and angular 
momenta of the mode equal to zero, as they should. Analogous mode in the n = 3 CSFCC 
cluster is depicted in panel (b) of Fig. El Again, the highest frequency mode is such that 
the central atom performs motion with the largest amplitude of all atoms in the cluster. 
The displacement patterns of the lowest frequency modes in n = 2 CSIC and n = 3 CSFCC 
clusters are depicted in panels (a) of Fig. |U and El respectively. It is of interest to note that 
the displacement pattern of the lowest frequency mode in n = 2 CSIC cluster is a sort of a 
"twisting" mode in which two halves of the cluster perform motions which look almost like 
the rotations around the same axis, but in opposite directions for the "upper" and "lower" 
halves of the cluster. 

Further inspection of the vibrational densities of states reveals significant differences 
between CSFCC and CSIC clusters even in the region of frequencies which contains the 
highest percentage of all vibrational modes (i.e. disregarding the high-frequency tail of the 
density of states in CSIC case). While p{oj) of CSFCC clusters obviously tends to the bulk 
(crystal) limit (except for the characteristic features around 10.8 ujq), the same does not 
hold for p(u) of CSIC clusters which behaves quite differently even for the largest cluster 
considered (n = 9). This was illustrated by a superposition of p(u) for CSIC and CSFCC 
n = 9 clusters in panel denoted by n = 9 in Fig. El Low frequency vibrations (up to about 
8 uj ) are quite similar in both types of sufficiently large clusters, but CSIC density of states 
does not exhibit characteristic features around 25 ujq, which in a FCC crystal (bottom panel) 
are a consequence of the van Hove type of singularity related to features of dispersion of 
longitudinal and transversal modes (phonons) at the crystal Brillouin zone edges (see e.g. 
Ref. E|). 
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The origin of the features around 10 u in CSIC and 10.6 ojq CSFCC clusters is in 
the large number of surface (poorly coordinated) atoms jl4j. Even for the largest cluster 
considered, made of 2869 atoms, 812 atoms are located at the cluster surface which makes 
about 28 percent (for smaller clusters the percentage is higher). This means that about 28 
percent of the vibrational density of states represents the modes in which surface atoms are 
significantly displaced from their equilibrium positions. Such modes are not present in the 
bulk crystal calculation. The peaks in p(uj) can in fact be related to the zone-edge frequencies 
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2611 . Additional 



of the Rayleigh-wave (RW) modes of surfaces of Lennard- Jones crystals 
confirmation of the identification of characteristic peaks can be found by inspecting panel 
(c) of Fig. 0] which displays the displacement pattern of a mode of n = 4 CSIC cluster 
with frequency of 10.33 u) . In this particular mode, mostly the surface atoms vibrate and 
the polarization vectors are oriented dominantly perpendicularly to the cluster surface. A 
similar pattern was found for larger CSIC and CSFCC clusters. However, the quality of 
the visual insight in the polarization pattern tends to degrade with the number of atoms 
in the cluster and this was the main reason for the choice of relatively small clusters for a 
visualization of the displacement patterns in Figs. |U and (see panel (c) in Fig. HJ). 
In Table ITT) the zero point energy, 

^ 37V-6 

f = -y: ku p (8) 
z P =i 

pertaining to the CSFCC and CSIC clusters is presented. These results are to some extent 
surprising. Although both the minimum and maximum frequency are smaller in the CSFCC 
clusters, the first moments of their frequency spectra are larger than in corresponding CSIC 
clusters, at least for the cluster shell numbers between 2 and 9 studied here. One can visually 
inspect that this is indeed so by looking again at the panel denoted by n = 9 in Fig. El This 
means that quantum ground state energies differ more than (classical) potential energies of 
the clusters. Thus, icosahedral closed shell clusters are even more preferred energetically 
when the quantum nature of the particles is important (in the cases when one can speak 
about the ordered ground state, when the approach presented here is adequate). To illustrate 
this effect and estimate its magnitude for clusters composed of rare gas atoms, in Fig. |U] I 
plotted the difference between the classical potential energies of CSIC and CSFCC clusters 
[AVp° = Vp (CSIC)— Vp (CSFCC), full circles] and the difference between the quantum ground 
state energies [AE = (CSIC) -V p ° (CSFCC) + F (CSIC)-F (CSFCC), empty diamonds] 
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TABLE II: Zero point energies (Fq) in units of hy / e/(Ma' 2 ) of Lennard-Jones icosahedral and 
cuboctahedral closed-shell clusters. The maximum cluster shell number is denoted by n. 
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12787.93 


13210.08 


6 


21138.61 


22359.53 


7 


33097.20 


34923.92 


8 


48960.54 


51457.37 


9 


68577.49 


72560.07 



as a function of cluster shell number n, and for clusters composed of Ne, Ar, and Xe atoms. 
These quantities cannot be obtained for general Lennard-Jones clusters in some reduced units 
of energy since the classical potential energy scales with e, while the zero-point energy scales 
with hujQ. The Lennard-Jonnes parameters used in this calculation were e(Ne) = 3.07 meV, 
a(Ne) = 2.75 A, e(Ar) = 10.35 meV, a(Ar) = 3.40 A, e(Xe) = 19.18 meV, a(Xe) = 4.1 
A 27]. This produces characteristic frequency scales fuj (Ne) = 0.29 meV, ?iu (Ar) = 0.31 



meV, and hu (Xe) = 0.19 meV. T 
as those obtained in Fig. 2 of Ref. 



'ie classical results (full circles) are practicaly the same 
It is obvious that quantum corrections become more 
important for lighter atoms. Even for Xe clusters, their effect is not negligible, especially 
when the differences between the CSIC and CSFCC clusters are considered. Intriguingly 
enough, the relative importance of quantum effects increases with the shell number. The 
behavior of free energy with temperature is such that the differences between the CSIC and 
CSFCC free energies increase in absolute value as the temperature increases. This means 
that at finite temperatures, the CSIC structure is even more favored than at T = K. The 
behavior of the sum of second two terms in Eq. (J1J) , 



3JV-6 * 3JV-6 

Fi = F-V p °= £ ^ + k B T £ In 
p=i z p=i 



cxp 



(9) 



k B T j 

with temperature is illustrated in Fig. [7|for n = 5 CSFCC (full line) and CSIC (dotted line) 
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clusters. Similar behavior is found for all n's examined in this work. 



IV. SUMMARY AND CONCLUSION 

The main result of this article is that quantal treatment of the low temperature properties 
of 2 to 9 shell CSIC and CSFCC Lennard- Jones clusters results in a larger binding energy 
of CSIC clusters. The difference between binding energies of CSIC and CSFCC clusters 
is larger than in classical treatment [^|. This results is somewhat surprising in the view of 
expectations put forth in the Introduction. It is plausible, although not shown in this article, 
that the presented quantal treatment would yield a quantum binding energy "crossover" 
between the CSFCC and CSIC clusters for larger shell numbers than obtained classically 
(14), i.e. that the CSFCC arrangement o: 
larger clusters than predicted previously 
the exact number of atoms for which the crossover takes place is dependent on the mass of 
atoms in the cluster. 



clusters would become energetically favorable for 
^]. On the basis of Fig. El it can be expected that 



[1] P. Moriarty, Rep. Prog. Phys. 64, 297 (2001). 

[2] J.J. Burton, J. Chem. Phys. 52, 345 (1970). 

[3] J.J. Burton, J. Chem. Phys. 57, 1980 (1972). 

[4] D.J. McGinty, J. Chem. Phys. 55, 580 (1971). 

[5] J. Xie, J.A. Northby, D.L. Freeman, and J.D. Doll, J. Chem. Phys. 91, 612 (1989). 

[6] J.A. Northby, J. Chem. Phys. 87, 6166 (1987). 

[7] B.W. van de Waal, J. Chem. Phys. 98, 4904 (1993). 

[8] J.P.K. Doye, D.J. Wales, and R.S. Berry, J. Chem. Phys. 103, 4234 (1995). 

[9] J.P.K. Doye and D.J. Wales, Chem. Phys. Lett. 247, 339 (1995). 

[10] B.W. van de Waal, Phys. Rev. Lett. 76, 1083 (1996). 

[11] R.H. Leary and J.P.K. Doye, Phys. Rev. E 60, R6320 (1999). 

[12] F. Calvo, J.P.K. Doye, and D.J. Wales, J. Chem. Phys. 114, 7312 (2001). 

[13] H.S. Lim, C.K. Ong, and F. Ercolesi, Surf. Sci. 269/270, 1109 (1992). 

[14] R. Meyer, L.J. Lewis, S. Prakash, and P. Entel, Phys. Rev. B 68, 104303 (2003). 



12 



[15] A. Siber, Phys. Rev. B 68, 033406 (2003). 

[16] W. Miehle, O. Kandler, T. Leisner, and O. Echt, J. Chem. Phys. 91, 5940 (1989). 

[17] S. Kakar, O. Bjorneholm, J. Weigelt, A.R.B. de Castro, L. Troger, R. Frahm, T. Mdller, A. 

Knop, and E. Riihl, Phys. Rev. Lett. 78, 1675 (1997) 
[18] K.R. Patton and M.R. Geller, Phys. Rev. B 67, 155418 (2003). 

[19] M. Born and K. Huang, Dynamical Theory of Crystal Lattices (Oxford University Press, 
Oxford, 1968). 

[20] H. Goldstein, Classical Mechanics (Addison- Wesley, 1980). 
[21] K. Huang, Statistical Mechanics (Wiley, New York, 1963). 

[22] D. Bressanini, M. Zavaglia, M. Mella, and G. Morosi, J. Chem. Phys. 112, 717 (2000). 

[23] R.E. Allen, G.P Alldredge, and F.W. de Wette, Phys. Rev. B 4, 1648 (1971); R.E. Allen, 

G.P Alldredge, and F.W. de Wette, Phys. Rev. B 4, 1661 (1971). 
[24] A. Siber, Phys. Rev. B 66, 235414 (2002). 

[25] N.P. Gupta and PK. Garg, Solid State Commun. 16, 607 (1975). 

[26] A. Siber, B. Gumhalter, A.P. Graham, and J.P. Toennies, Phys. Rev. B 63, 115411 (2001). 
[27] G. Stan, M.J. Bojan, S. Curtatolo, S.M. Gatica, and M.W. Cole, Phys. Rev. B 62, 2173 
(2000). 

FIG. 1: Closed shell icosahedral (CSIC, top panel) and cuboctahedral (CSFCC, bottom panel) 
Lennard-Jones clusters in their minimum potential energy configuration for maximum shell num- 
bers n =2,3, and 4 (from left to right). 

FIG. 2: Vibrational densities of states of CSIC clusters as a function of maximum shell number, 
n = 2,. ..,9. The vibrational density of states of Lennard-Jones FCC crystal is displayed in the 
bottom panel for comparison. 
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FIG. 3: Vibrational densities of states of CSFCC clusters as a function of maximum shell number, 
n = 2,. ..,9. The vibrational density of states of Lennard-Jones FCC crystal is displayed in the 
bottom panel for comparison. In panel denoted by n = 9, the vibrational density of n = 9 CSIC 
cluster is also displayed and denoted by a dotted line. 

FIG. 4: Eigenvectors (displacement patterns) of three vibrational modes pertaining to CSIC clus- 
ters. The modes depicted in panels (a), (b), and (c) are those denoted by a, b, and c in Fig. [21 
respectively. Equilibrium positions of cluster atoms are denoted by small cubes. The displacement 
vectors are multiplied by 10 in panels (a) and (b), and by 20 in panel (c). 

FIG. 6: The difference between the classical (full circles) and quantum (empty diamonds) ground 
state energies of CSIC and CSFCC clusters as a function of maximal shell number for Ne (top 
panel), Ar (middle panel), and Xe (bottom panel) clusters. 



FIG. 7: Free energies of n = 5 CSFCC (full line) and CSIC (dotted line) clusters measured from the 
corresponding classical potential energy minima [see Eq. Q] as a function of reduced temperature 
k B T. 



FIG. 5: Eigenvectors (displacement patterns) of some of the vibrational modes pertaining to CS- 
FCC clusters. The modes depicted in panels (a), (b), and (c) are those denoted by a, b, and c 
in Fig. EH respectively. Equilibrium positions of cluster atoms are denoted by small cubes. The 
displacement vectors are multiplied by 15 in panels (a),(b), (c). 
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